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The phase diagrams and superconducting properties of the extended Hubbard model with pair 
hopping interaction, i.e. the Penson-Kolb-Hubbard model are studied. The analysis of the model 
is performed for d- dimensional hypercubic lattices, including d — 1 and d — oo, by means of the 
(broken symmetry) Hartree-Fock approximations and, for d — oo, by the slave-boson mean- field 
method. For d = 1, at half-filling the phase diagram is shown to consist of nine different phases 
including two superconducting states with center-of-mass momentum q = and q = Q (rj- pairing), 
site and bond-located antiferromagnetic and charge-density wave states as well as three mixed phases 
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with coexisting site and bond orderings. The stability range of the bond-type orderings is shrank 
with increasing lattice dimensionality d and for d = oo the corresponding diagram consists of four 
^ ■ phases only, involving exclusively site-located orderings. Comparing the pair hopping model with 

the attractive Hubbard model we found in the both cases gradual evolution from the BCS-like limit 
to the tightly bound pairs regime and a monotonic increase of the gap in the excitation spectrum 
with increasing coupling. However, the dynamics of electron pairs in both models is qualitatively 
different, which results in different dependences of condensation energies and critical temperatures 
on interaction parameters as well as in different electrodynamic properties, especially in a strong 
— ■ coupling regime. 
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The purpose of the present work is the analysis of phase diagrams, electronic orderings and superconducting 
properties of the extended Hubbard model with pair hopping interaction, i.e. the so-called Penson-Kolb-Hubbard 
55 ■ (PKH) model, 

- J Ei,/4r c a c J\i c 3T - mE*, ct (l) 



'■^J ■ where the prime over the sum means restriction to nearest neighbor (n.n) sites, t denotes the single electron hopping 
integral, U is the onsite density-density interaction, J is the pair hopping (intersite charge exchange) interaction and 
/i is the chemical potential. In the absence of the U term the Hamiltonian (Q) reduces to the Penson-Kolb (PK) 
H ; modelJa 

£> . We will treat the parameters t, U, J as the effective (phenomenological) ones, assuming that they include all the 
possible contributions and renormalizations like those coming from the strong electron-phoiion couplings or from the 
coupling between electrons and other electronic subsystems in solid or chemical complexes^ (such that the values of 
U and J can be effectively either positive or negative). It is notable that formally J is one of the off-diagonal terms 
of the Coulomb interaction — J = (ii\e 2 /r\jj)fl describing a part of the so-called bond-charge interaction, and the 
sign of the Coulomb-driven charge exchange is te]?ically negative(repulsive, J < 0). However, the effective attractive 
interaction of this form (J > 0) is also possibldju and in particular it can originate from the coupling of electrons 
with intersite (intermolecular) vibrations via modulation of the hopping integral^ or from the on-site hybridization 
term in a generalized periodic Anderson modeLuu 

The PKH model is one of the conceptually simplest phenomenological models for studying correlations and for 
description of superconductivity of the narrow band systems with short-range, almost unretarded pairing. It includes 
a nonlocal pairing mechanism (the pair hopping term J) that is distinct from the on-site interaction in the attractive 
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Hubbard model and that is the driving force of pair formation and also of their condensation. Thus, the supercon- 
ducting properties and the evolution from the Cooper pair regime to the strong coupling local pair regime can be 
essentially different in these two models. 

While most of the basic properties of the attractive Hubbard model seems to be at present well unders 
several years of intense studies, the PKH model has been investigated only in a few particukuL-Jimitsl 
main efforts concerned the ground state phase diagram of the half-filled one dimensional PKHEllj'tEl and PKL 
models. In the case of the PKH model these problems were studied by both, momentum-space renprmalization-group 
(MSRG) and the finite-size (exact diagonalization of finite-size cells) methods (for U,J> O^ElJay the real space 
renormalization-group (RSRG) (for U > 0)]B, by the continuum-limit field theory (CFT) approacblij (for U > 0) and 
within the Green's function formalism in the mean-field approximations. However, in all these studies, except [ [l5| , 
the possibility for the bond- located orderings was not considered and the exact form of the phase diagram in the whole 
range of parameters — oo < U/t, J/t < oo has not been established up to now. The properties of the PKH model for 
higher dimensional lattices (1 < d < oo) and aphitrary electron concentration (0 < n < 2) have not been studied yet, 
except for the limiting case of zero bandwidth.t3 The latter limit was analyzed by the variational approach, in which 
the U term is treated exactly and the intersite J term - within mean-field approximation!!^ (such an approach yields 
exact results for d = oo). 

In the paper we will study the PKH model for the case of d-dimensional hypercubic lattices (1 < d < oo) and 
arbitrary, positive as well as negative, U and J. In the analysis we will apply a broken symmetry Hartree-Fock 
approximation (HFA) (Sec. 2) supplemented, for d — oo, bv-the slave boson mean-field approach (SBMFA) (Sec. 3). 
In the case of the Hubbard model and its various extensionscl the former approach is known to give credible results at 
T = for any U as far as the energy of the ground state and energy gap in the ordered states is concerned. It usually 
provides qualitatively correct ground state phase diagrams for arbitrary dimensions if all the proper broken symmetry 
phases are included into the analysis. Moreover, for the electronic models with intersite interactions only, the HFA 
becomes an exact theory in the limit of infinite dimension (d = oo). At T > the HFA is much less reliable, especially 
for low dimensional systems and the limits of strong coupling, as it neglects short-range correlations and the effects 
of collective excitations. An obvious weakness of the HFA (both at T = and T > 0) is inadequate description of the 
normal (nonordered) phase. This failure is a consequence of the fact that the HFA greatly overestimates the energy 
of the phases without long-range order. Going beyond the HFA we will use the SBMFA. The slave-boson method 
is in principle not restricted to weak ox-.strong coupling and it is an improvement over the former treatment since 
it takes into account local correlations .ED We will apply the SBMFA only for d — oo, where the intersite coupling J 
can be treated adequately. For finite dimension (d < oo) the SBMFA treatment of intersite interactions is technically 
involved and to our knowledge it has not been analyzed consistently so far. 



II. GENERAL FORMULATION AND THE HARTREE-FOCK ANALYSIS 



In the system considered several types of superconducting, magnetic and charge orderings can develop. In 
the following we will study the case of alternating (hypercubic) lattices with nearest-neighbor single electron 
hopping t and pair hopping J, and restrict our considerations to the one- and two-sublattice orderingsJ13 de- 
scribed by the following order parameters: the superconducting with the s-type (S) and the retype (77) pair- 
ing: x s = ^Ei( c a c *T> = W Efc( c -fc! c fcT)> x v = ^E l elQR, ( c a c *T> = ^E/c( c -fc+Qi c fcT>; the antiferromag- 
netic (AF) with the staggered magnetization located on sites (sAF) or on bonds between sites (bAF): x s af = 

^£i,** eiQR< ( c L c i*>» x bAF = ^J2'i,j>i,a^ eiQRi (4a c ^) = mT,k,a a Vk(ct a c k +Qa); the charge density wave 
(CDW) with the on-site (sCDW) or the bond zigzag (bCDW) modulation of charges: x sC dw = jjy J2i a e lCJRi (c\ a c icr ), 

x b cDW = 27v T!i,j>i, a e ' iQR * ( c L c jv) = ITv Efe, C r 7 7fe< c L c fc+Q^) J where Vk = iJ2a sin ( ka ) and Q = (n/a, n/a, ...). We 
assume that the sites are ordered in an ascending way along the crystallographic axis and for the case of the bond 
zigzag parameters the sum is restricted to the nearest neighbor sites j, which followed the i-th site. The number of 
electrons per lattice site is given by n = j^^i (r^icr ^) ■ ^ n * ne case °f the AF phase we quoted above only sAF 2 
and bAF 2 orderings, corresponding to a z-component magnetization located on sites and bonds, respectively, and we 
omitted s(b)AF x , s(b)AF !/ . Due to the SU(2) - spin symmetry of the PKH model the latter orderings are strictly 
degenerated with s(b)AF 2 . 

Within the framework of the broken-symmetry Hartree-Fock approach the mean-field Hamiltonian in the momentum 
space fc, including all types of orderings is given by 

Hhf = ^{tk - M + \ n ) c \o c k<? + (U - zJ) ^2(x s cl^_ kl + h.c.) 

k,(T k 
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-(U + zj) J2^ C U C -k+Ql + h - C -) ~ Ux sAF aC ka C k+Q<J 
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XbCDW 53 % c L c fc+Q<r > ( 2 ) 



where t k = —t^k, t = t + 2pJ/z, z is the number of nearest neighbor sites (for the hypercubic lattice of c?-dimension: 
z = 2d), and p denotes the Fock term: p = ± T,'i,j,A c L c j^) = IF £fe, CT 7fc(4<r c ^)> with 7fe = E Q cos(fca). 
The eigensolutions of the Hamiltonian (0) and the corresponding free energy 



F = ~ ln[Tr{exp(-^ F )}] + (H - H„ F ) HF + f iN e , (3) 

where A-t= l/fc^T and iV e denotes the number of electrons in the system, can be determined by the standard 
methodsE^I with either the Green's function or the equation of motion approach. If the solutions corresponding to the 
pure phases (i.e. the phases with only one type of order) are analyzed, the free energy (||) may be expressed in terms 
of the eigenvalues of Hhf in the form 

— =-p(n-l) + —n 2 + -Jp 2 +A a \x a \ 2 
N 4 z 



L^l n [2cosh(^)], (4) 



where r = ±, Jt = /i — Un/2, V a is an effective coupling strength for the a phase, which is Vs = —U + zJ, 
V n = —U-zJ, V s af = U, V b AF = —2J/z, V sC dw — —U and V bC DW — 2J/z , A a — V a for a — S, 77, sAF, sCDW and 

A a = 2V a for a = bAF, bCDW. The electronic spectrum is Ef k = ± v / (e k - ~p) 2 + Vgx' z s , E± k = e k ± ^//J 2 + V 2 x^, 

E fAF k = M ± + VsAf x 'Iaf> 4f k = 7* ± V e fc + V h 2 AF \r, k x bA F\ 2 , Ef CDW k = ~p ± V e fe + V sCDW x icDW and 
E b ± CDW k =~p± \Zel + V b 2 CDW \r]kXbCDw\ 2 for the S-, the 77-, the sAF-, the bAF-, the sCDW- and the bCDW phases, 
respectively. In the derivation of the eigensolutions we have assumed an alternated lattice, i.e. £fc+q = — (-k- 

For arbitrary electron concentration n the stable solutions are determined as the minimum of F with respect to the 
variational parameters x a (a = S, 77, sAF, bAF, sCDW, bCDW), p and p, i.e. by the equations 

dF/dx a = , dF/dp = , dF/dp = . (5) 

Besides the pure phases there are also solutions for various mixed type orderings. We have analyzed the stability 
conditions for all such states and found that some of them can be stable in a definite range of parameters. They 
are summarized in the Table I together with the corresponding order parameters. For example, we present here the 
equations describing the mixed s+bAF phase. In this case the free energy ^ is expressed in terms of the eigenstates 
as 

^ = Ji(n - 1) + ^-n 2 + -Jp 2 + V s afx 2 s af 
N 4 z 

+2V hAF \x hAF \ 2 -^5>[2cosh( ^^* )] (6) 



fc,-, 



with the electronic spectrum Ef bAF k = 7J± ^/e k + V 2 AF x 2 AF + V b 2 AF \r)kXbAF\ 2 and x s af, x b AF, P and \x are deter- 
mined by a set of self-consistent equations: dF/dx s AF — 0, dF/dxbAF — 0, dF/dp = and dFjd\x = 0. 

In order to determine the mutual stability of the phases considered one has to find all the possible solutions 
and compare the corresponding free energies. In the weak and strong coupling regimes we were able to derive several 
analytical expressions concerning the energy gaps, the order parameters and the critical temperatures, but in a general 
case numerical methods had to be used. At T = we performed complete numerical analysis of all the solutions for 
the whole range of the parameter values and the resulting phase diagrams are presented in Fig.|l| for the ID chain and 
the hypercubic lattice of the dimension d — 00. The renormalized parameters are: J* = Jd and t* = t^fd. 

For d — 00 the density of states (DOS) is p(e) = exp[— e 2 /(8t* 2 )]/(v / 87ri*). In this case there are not stable states 
with the bond type of ordering as all bond parameters disappear in the limit d — -> 00. Also, the Fock term p is then 
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irrelevant as the effective width of the electronic band W e ff = 4fd = At*^/d + 4J*p/d and the second term disappears 
for d = oo. This is in contrast to the d = 1 case (Fig.|j]a), where the AF and CDW orderings of the bond type can 
exist in a wide range of parameters (the former for J < and the latter for J > 0). The bond type ordering can 
also coexist with the on-site type ordering, as it is seen in Fig. pi for the mixed s+bAF phase. There have been also 
found very narrow regions of the stable mixed phases: bAF+sCDW (for U < 0) and sAF+bCDW (for U > 0). The 
curves separating the sAF- and the bAF-type orderings are the lines of second order phase transition, at which the 
parameter x s af or XbAF disappears. In the lattices of dimension 1 < d < oo one can analyze more complex bond 
orderings (e.g. the phase of fluxes), however, the ranges of stability of all the bond-ordered phases will be gradually 
shrank with increasing lattice dimension. 

The J* dependence of the order parameters for U/t* = 3 is presented in FigJ|, where the upper part is for the 
d = 1 system and the lower part for d = oo. Fig.|2|a shows a wide range of the mixed s+bAF phase with x s af ^ and 
XbAF ^ 0. The parameter XbAF —> for J* /t* — > —2.83, indicating the second-order transition. In the case presented 
in Fig.|a (U/t* = 3) the mixed sAF+bCDW phase is stable only in a very narrow range 2.09105 < J* /t* < 2.09425. 



III. SLAVE-BOSON STUDIES 

In the previous paper0 we showed that the slave boson mean-field approach (SBMFA) gives reliable results for the 
ground state properties of the attractive Hubbard model in the whole range of coupling \U\ and arbitrary electron 
concentration n. Therefore, we also applied this method to the present model (|l|). As the SBMFA takes into account 
the onsite electron correlations and neglects the short-range intersite correlations (the Fock term and the bond type 
orderings are omitted), we have concentrated on the case of d = oo lattice, where the mean field treatment of intersite 
interactions becomes exact. 

In the slave-boson approach each local state is described by a fermi operator fi a and two types of bose operators Pi 
and bi, which correspond to two vector fields: a field of local magnetic moments and that of local charges. The com- 
pleteness condition means that length and direction of the vectors pi and bi can vary from site to site, but a sum of . thei r 
length is always pf + b 2 = 1. We use the spin- and the charge-rotationally invariant slave-boson representationEJIIZl, 
in which the order parameters are expressed by xs = jfY.i. P ( b L^ + b l y b iy)i x v = w Xrf.p e jQR * (b\ x b ix + b\ y bi y ), 

XsCDW = jf J2i e lQRi (bl z b iz ) , and x sA f = j? J2i e ' QRi (PiJPiz) > for tne superconducting, sCDW and sAF phase, 
respectively. In the mean field studies we confine ourselves to the temperature T = and neglect space and time 
fluctuations of the bose fields. The operators pi and bi are replaced by their expectation values, which in the following 
are treated as variational parameters. The SBMFA is, therefore, a variational method on a trial state described by 
the Hartree-Fock wave functions (being equivalent to the Gutzwiller approximation). Ell The free energy is the sum of 
the fermionic and bosonic parts, and for the a phase, where a = S, 77, sCDW, sAF, and any given n it can be written 
in the following unified form: 

pSBMFA p p b 

N = ~N + ~N ~ 

k,r 

+C a - (\ + + 25) - 2\ a x a , (7) 

where r = ±, C s = -2J*x%, C n = 2J*x%, C sCDW = 0, C sAF = 0, b 2 = {b\ x b lx + b\ y b ly + bjb lz ) and 25 = 
n — 1 = (b\ z b iz ), Ao and X a are the Lagrange multipliers. The fermionic spectrum is Ef^^ = ±y/ (qs^k + Ao) 2 + As, 

E? k ± = q n e k ± + Ef c % w k = -Ao ± y/qsCDWet + KcDW and E sAF k = ~ A ± VlsAF4 + X'i AF - Its 

fc-dependence is analogous to that obtained in the HFA with the bandwidth reduced by the factor 



_ 2p\tf + ^ - 4x1 - 4^ 2 ) 



for a — S, 77, sCDW, and 



2bV + ^_ Axj AF -A5 2 ) 
qsAF = 1-4x^-4* • (9) 



The stable solutions are determined from the minimum of the free energy p SBMFA with respect to x a , A Q , Aq and b. 
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In determination of the phase diagram of the half-filled PKH model we compare the free energies J 
responding to the S, 77, sCDW and sAF phases. Their values are different from those obtained in the HFA and 
depend on the band narrowing factors q a . These factors are important parameters. In the ipjmaLphase [x a = 0) 
the band narrowing process can lead to insulating phase (qw = 0) for large coupling \U\ ^> i.E30'EJE3 However, for 
alternating lattices and n = 1 the normal phase is not a ground state as its free energy is always higher than that 
of the long-range ordered phases (S, 77, sCDW and sAF). For all these phases the band narrowing factors q a are 
close to unity, for example, in the attractive Hubbard model: 0.954 <(7s<lin<i = oo( see also Ref.[ p2[). Thus, 
the SBMFA free energies of the ordered phases are relatively close to the corresponding HFA results. The SBMFA 
phase diagram of the PKH model for d = 00 is, therefore, very similar to that given in Fig.^3. In particular, the 
location of the S-sAF phase boundary in the ground state phase diagram can be expressed most conveniently in terms 
of the deviation e c = J* /U — 1 from the line J* /U = 1. Within the HFA, the fhsAF phase boundary is given by 
J*/U = 1 for any t* , and for d = 00 it agrees with a rigorous solution at t* — O.ta Within the SBMFA, e c is found 
to depend sensitively on the strength of the interactions and one obtains that e c > for any 00 > t* > 0, with a 
maximum deviation e c « 0.02 found for U/t* w 3.5 and with e c — ► for t* — > as well as t* — > 00. It means that the 
hopping term slightly extends the stability range of the sAF phase with respect to the S phase. Notice that similar- 
results are obtained for the extended Hubbard model with nearest neighbor density-density repulsion W*. In that 
case Monte-Carlo simulationsES and perturbational treatmentsE.il show that for t ^ the actual phase boundary is 
also slightly shifted upward relative to the line W* /U — 1 predicted by the HFA. 

Although the SBMFA gives minor changes in the ground state energies, other physical characteristics are modified 
in a much more pronounced way. We will show it analyzing the gap E g in the excitation spectrum determined within 
the SBMFA as well as the HFA. Fig.|| shows dependences of E g on J* in the case of U/t* — 3 and n = 1. The value 
of Eg BMFA is reduced with respect to E^ FA . The results are closer to each other for larger couplings |J*|, where 
the onsite correlations become less relevant. The maximum reduction is seen for the gap in the r\ state, which at the 
transition line is reduced by a factor j n = E g ® MFA /E g HFA = 0.56. Fig.^ presents the U/t* dependence of E g for 
J* ft* — —2 and n = 1. The gaps E g s af and E g s cdw corresponding to the sAF and sCDW phases do not depend 
on J*, and they are the same as in the usual Hubbard model (J* =0). In the lower part of Fig.^ the reduction 
parameter 7 Q is shown. The minimum value of -f v is 0.11 for E g ^ MFA close to the transition point to the sCDW 
phase. The energy gaps E^^ IFA and E^^ 1 ^ are maximally reduced for a weak coupling |L//£*| <C 1. In this limit 
one can find that 



and 



E™ A =At*e XP [-^f] ,10. 



j^SBMFA o 

7a = g ° =cxp[--ft = 0.554855, (11) 

where a = sAF, sCDW and S, A = 4 v / 2V' y / 2 = 4.23871 and 7 = 0.577216 is Euler gamma constant. (The value @ 
is larger than 75 = exp [—3/4] = 0.472267 obtained for the rectangular density of states.HH ) 



IV. DISCUSSION AND CONCLUDING REMARKS 



Let us compare the properties of superconducting phases and their evolution with a change of coupling and con- 
centration for the three limiting cases of the model (|l|): i) the attractive Hubbard model with U < 0, J — , ii) the 
PK model with U = 0, J > 0, and iii) the PK model with U = 0, J < 0. We will discuss qualitative differences and 
similarities in the behavior of the system for these limits and stress distinct features of each case. 

In the first two cases the pairing interaction favors the on-site s-wave superconductivity (S), whereas in the third 
one - the 77-pairing. Moreover, the later two cases include a nonlocal pairing mechanism (J) that is distinct from the 
zero-range instantaneous interaction existing in the i) case. The difference between i) and ii) occurs in the case of the 
half-filled band. At n = 1 the U < Hubbard model posses SU(2) symmetry of the charge sector and is characterized 
by coexistence of the sCDW and the S ordering (these phases are strictly degenerated) in the ground state. No such 
degeneracy occurs in the PK model, as its charge sector is governed by the U(l) symmetry for any n. 

For i) and ii) at T = the S-phase is stable for any nonzero interaction (U < or J > 0) and arbitrary n 
(0 < n < 2). In both these cases the evolution of the S-phase from the BCS like superconductivity, with extended 
Cooper pairs, to superconductivity of composite bosons (local pairs) with increasing coupling is continuous. At T = 
the appropriate boundary between both regimes can be located (after Leggett [ |5|) from the requirement that the 
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chemical potential in the superconducting phase reaches the bottom of the electronic band, i.e. from = — W/2. 
For d = 2 lattice the borderlines as a function of n are shown in Fig |. As we see for both models with increasing 
n the boundaries are shifted towards higher values of coupling. For d = 3 the corresponding plot has qualitatively 
similar form, except n — > limit, where there is a critical value of coupling for pair formation. 

For the case iii) the 77-phase is stable only below a critical value of J and for d < 00 the local pair regime is 
reached directly after crossing the r\ phase boundary. The critical value J c depends on the lattice structure, the 
lattice dimensionality (d) and the band filling (n). The estimations of J c for various cases are collected in Table II. 
Except d — 00, the transition at J c is of the first order and characterized by an abrupt change in the structure of the 
ground state. For d = 00 the phase stable for J c < J < is a normal metal without any long-range ordering (for 
any n), whereas for d < 00 and n — 1 that phase is insulating and antiferromagnetic with bond- type modulation of 
magnetization. 

The evolution of the gap parameters x a (a = S, r\ and bAF) with increasing interaction V a (for all three cases) 
is presented in Fig.^J for d = 1 (Fig.^a) and d = 00 (Fig.^Jb). The corresponding plots for d = 2 and d — 3 lattices 
have qualitatively the same form to those for d = 1. For the sake of comparison we have shown also the SBMFA 
results (curves with diamonds in Fig.||) calculated for the attractive Hubbard model in d = 1 and d — 00. Notice 
the first order transition from the bAF state to the 77 state in the J < PK model for d = 1. In this case x v has 
the maximum value 1/2. On the contrary for d = 00 the phase transition to the 77-state is of a second-order, and x v 
continuously increases with decreasing J* (it never saturates for a finite J*). For the J > PK model and the U < 
Hubbard model one observes a continuous evolution of the order parameter x$ with increasing V a and an exponential 
(BCS-like) behavior of xs in the weak coupling limit. 

As we have already pointed out, for both models there is a crossover from the BCS-like limit to the tightly bound 
pairs regime with increasing coupling and this evolution of the superconducting (S) phase is gradual. However, 
the thcmiodynamic and electromagnetic properties of both the models are very different beyond the weak coupling 
limit.E3c£l To illustrate the situation we have plotted in Figj^ the condensation energies, i.e. the difference of the free 
energy in the normal and in the superconducting phase, AF = F N — F s , as a function of the coupling parameters 
V a . These results have been obtained for the d = 1 chain at n = 0.8, but for the lattices of other dimension and 
other electron concentrations one gets qualitatively similar dependences. They are in good qualitative agreement with 
results of perturbational expansions Jor-the models considered both in the weak coupling (V a /t <C 1) as well as in the 
strong coupling regimes (V a /t ^> l).EZlc3 As the square of the thermodynamic critical field H 2 is proportional to AF, 
we conclude that in the attractive Hubbard model this quantity (similarly as the critical temperature T c ) increases 
exponentially for small values of \U\, then it goes through a round maximum and it decreases as t 2 /\U\ for large \U\. 
On the contrary, in the PK model we found no maximum of H 2 and T c at intermediate coupling and both these 
quantities increase linearly with J for large J. Also the behavior of the penetration depth Al and the pair mobility 
t p is different. In the strong coupling limit \ 2 L oc l/t p increases with \U\ in the former model (A^ oc \U\/t 2 ), while 
it decreases with J in the PK model (\ 2 L cx l/J).pEjb is in agreement with studies of collective excitations performed 
using a generalized random-phase approximation.!!-] The collective-mode velocity increases with J in the PK model, 
in contrast to the attractive Hubbard model where it decreases with the coupling \U\. 

The phase diagram of the half-filled one-dimensional PK model (U = 0, J 7^ JO) derived within the HFA is in 
agreement with that obtained by the density matrix renormalization group methodJ3 : E3il3 For J > both approaches 
predict a continuous second-order transition to usuaLs-wave pairing state at J = + , with no additional transition 
for any J > (in contrast to the earlier predictionsErliil) . We have found that (at least for alternating lattices), this 
phenomenon remains unchanged in higher dimensions (including the exactly solvable case of d = 00) and does not 
depend on the band filling. With increasing coupling there is a gradual crossover from the BCS-like superconductivity 
to the superfluidity of tightly bound local pairs. On the contrary, for J < the HFA predicts that the 7y-phase is 
stable only above a critical value of \ J\ and that the transition at J c is of the first order for any d < 00. Let us stress 
that for n = 1 the values of J c calculated within the HFA ape, p.A£ery good quantitative agreement with the results of 
other more elaborated treatments available for d — 1 chairJi3'tlrll3(cf. Table II). Moreover for n — > the HFA yields 
exact results for J c for any dimension. 

We have found that the interplay between the on-site Coulomb interaction U and the intersite pair hopping J 
in the PKH model can stabilize several new ordered phases absent in the usual Hubbard model (U ^ 0, J = 0) 
and in the usual PK model (U = 0, J ^ 0): the bCDW and the mixed bCDW+sAF phases (for U > 0, J > 0), 
the mixed s+bAF phase (for U > 0, J < 0) as well as the sCDW and the mixed bAF+sCDW phases (for U < 0, 
J < 0). The new phases predicted by our broken symmetry HFA approach (which can be truly long-range in d > 2) 
indeed need further examination by more rigorous methods as the exact diagonalization of small systems, density 
renormalization group, etc. We should point out however that our findings concerning the bond-ordered solutions 
are clearly supported by recent work of Japaridze and Miiller-HartmanrtJ performed for the d = 1 PKH model with 
U > in the weak coupling, using the continuum limit field theory approach and bosonisation technique. In contrast 
to previous studies ,Bl3'tll which have not considered the possibility of bond- located orderings, present results and 
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those of Ref. [ |l5| indicate that the bCDW (bAF) state but not the sAF or sCDW phases, is unstable with respect to 
transition into the S (tj) phase with increasing J (—J) for U > 0. 

We have also compared the superconducting properties of the PK model with those of the attractive Hubbard 
model. Although the energy gaps have similar dependences on the coupling parameters in the both models (see also 
[ ^)), dynamics of electron pairs is qualitatively different, which results in different electrodynamic properties and 
different coupling dependences of T c , especially in a strong coupling regime. 
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TABLE I. Phases considered and the corresponding order parameters. 



Type of phase 



n 

sAF 
bAF 
s+bAF 
sCDW 
bCDW 
s+bCDW 
bAF+sCDW 
sAF+bCDW 
S+sCDW 
S+bCDW 
ry+sCDW 
77+bAF 
?7+sAF+bAF 



Order parameters 



X S AF 
XbAF / 
XsAF =fi 0, X bA F =fi 
XsCDW / 
XbC'DW 7^ 
XsCDW 7^ 0, IfcCDH' 
XbAF 7^ 0, XsCDIV 7^ 
IsAF / 0, XbCDW # 

a;s / 0, Iscdw # 

IS / 0, ItCDlC / 

/ 0, Ijcaw / 

X v / XbAF 
I,, 0, X S AF 7^ 0, XfcAF / 



TABLE II. The HFA estimates of the critical value of J below which the 77-state has lower energy than the normal state 
in the PK model with J < 0. In the limit n — > the exact solution is J* =r-rW — — 2\[dt* . The results obtained b% 
density matrix renormalization groupEj (DMRG), Lanczos diagonalizatiorcj and real space renormalization-groupliil 
(RSRG) methods for the ID chain are given as well. 



system 


rjr 


n = 1 


n — > (exact) 


d = 1 


-tt/2 « -1.5708 
-1.5 [0 , -1.75 [|ll], -1.65 [0 


-2 


d = 2 square lattice 
rectangular DOS 




-8/3^- -1.8856 


-2-Z2 « -2.8284 


d = 3 sc lattice 
elliptic DOS 


-1.7028 
-3v^tt/8 « -2.0405 


-2^ ~ -3.4641 


d = 00 


-v 7 ^- -2.5066 


—00 
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FIG. 1. Phase diagram of the half-filled PKH model for the ID chain (Fig. a) and for the d = cx)-hypercubic lattice (Fig.b) 
determined within the broken symmetry HFA. The region of the mixed AF state is denoted by sbAF. Close to the boundary 
lines separating the sAF and bCDW states as well as the sCDW and bAF states there are very narrow regions (narrower than 
thickness of the curves in the figure) of the stable mixed ordered phases (sAF+bCDW, for U > and sCDW+bAF for U < 0). 
First-order and second-order transition phase boundaries are marked by solid and dashed curves, respectively. The SBMFA 
phase diagram for d = oo is almost identical with Fig.b (see discussion in Sec. 3). 
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FIG. 2. The J* /t* dependence of the order parameters in the ground state of the d = 1 (Fig. a) and d — oo (Fig.b) system, 
for n — 1 and t//i* = 3. The stability ranges of the different phases are indicated by the vertical dashed lines. The mixed 
sAF+bCDW state (Fig.a) exists for J* /t* £ (2.09105,2.09425). 
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FIG. 3. The dependence of the gap in the excitation spectrum on J* /t* calculated in the SBMFA (solid curves) and in the 
HFA (dashed curves) for d = oo hypercubic lattice, U/t* = 3 and n = 1. The stability ranges of the different phases are 
indicated by the vertical dashed lines. 
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FIG. 5. Boundary lines between the regions of the BCS-like and the local pair superconductivity calculated as a function of 
n in the case of d = 2 lattice for the U < Hubbard model (V a = —U) - solid curve, as well as for the PK model with J > 
(V a = 4J ) - dashed curve, and with J < (V a = — 4J) - dotted curve. For the PK model with J < the »?-type Cooper pairs 
are stable only above the long-short dashed. 
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FIG. 6. The plots of the gap parameters x a = E™ /V a as a function of the coupling parameters V a for the S phase of the 
attractive Hubbard model (dashed-dot curve, V a = —U) and the PK model (solid curve, V a = 2J*) as well as for the r\ phase 
(dotted curve) and the bAF phase (dashed curve) of the PK model (V a = —2 J*), calculated within the HFA, in the case of 
d — 1 (Fig. a) and d — oo (Fig.b) lattices, n = 1. For d = co the plots of xs vs. V a for the PK model and the U < model 
have the same form (solid curve). For the sake of comparison the SBMFA results for the U < Hubbard model in d = 1 and 
d = oo are presented by the curves with diamonds. 
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FIG. 7. Difference between the free energies for the normal and the superconducting phase plotted as a function of coupling 
parameters V a for the attractive Hubbard model (dashed-dot curve, V a = —U) and the PK model (solid curve, V a = 2J) as 
well as the difference of the free energies for the normal and the ?7-type superconducting phase for the PK model vs. V a = —2 J 
(dotted curve). The derivations were performed for the d = 1 chain at n — 0.8 and T = by means of the SBMFA for the 
attractive Hubbard model and by the HFA for the PK model. 
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